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THE SPIN- TEMPERATURE THEORY OF DYNAMIC NUCLEAR 
POLARIZATION AND NUCLEAR SPIN-LATTICE RELAXATION 1 

By Charles E. Byvik and David S. Wollan^ 

Langley Research Center 

SUMMARY 

A detailed derivation of the equations governing dynamic nuclear polarization (DNP) 
and nuclear spin- lattice relaxation by use of the spin- temperature theory has been carried 
to second order in a perturbation expansion of the density matrix. Nuclear spin diffusion 
in the rapid diffusion limit and the effects of the coupling of the electron dipole-dipole 
reservoir (EDDR) with the nuclear spins are incorporated. The complete expression for 
the dynamic nuclear polarization has been derived and then examined in detail for the limit 
of well-resolved solid-effect transitions. Exactly at the solid-effect transition peaks, the 
conventional solid-effect DNP results are obtained, but with EDDR effects on the nuclear 
relaxation and DNP leakage factor included. Explicit EDDR contributions to DNP are dis- 
cussed, and a new DNP effect is predicted. 

INTRODUCTION 

Soon after the classic work of Bloembergen, Purcell, and Pound (ref. 1), Purcell 
(quoted by Bloembergen in ref. 2) suggested that the spin temperature of paramagnetic 
ions in a solid could be raised above the temperature of the lattice by a strong radio- 
frequency field at the resonance frequency of the electronic spins and, furthermore, that 
the neighboring nuclear spins would follow the electronic spin temperature. This process 
of driving the nuclear spin temperature away from the temperature of the lattice by pump- 
ing the electron spins suggested the possibility of dynamic nuclear polarization (DNP). 

In 1953, Overhauser (refs. 3 and 4) demonstrated theoretically that the spin temper- 
ature of nuclei in a metallic solid could be lowered through the contact hyperfine interac- 
tion of the nuclei with the conduction electrons when the electron spin resonance transition 

Ipart of the information presented herein was included in a thesis entitled 
"Dynamic Nuclear Polarization in Samarium- Doped Lanthanum Magnesium Nitrate" sub- 
mitted by Charles E. Byvik in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in Physics, Virginia Polytechnic Institute and State University, 
Blacksburg, Va., September 1971. 

^Assistant Professor of Physics, Virginia Polytechnic Institute and State University, 
Blacksburg, Va. (Consultant at NASA Langley Research Center). 



is saturated. The experimental verification of Overhauser's predictions of DNP in met- 
als was reported in 1956 (ref. 5). The idea of dynamic polarization of nuclei was soon 
extended to nonmetallic materials (refs. 6 to 9). The theory of DNP and nuclear spin- 
lattice relaxation in ionic crystals with paramagnetic impurities is the subject of the pres- 
ent investigation. 

Two theoretical approaches have been used to describe the spin dynamics of such 
systems: the rate-equation approach (refs. 10 to 14) and the spin-temperature theory 
(refs. 15 to 23). The rate-equation approach assumes that the paramagnetic electron spins 
and the nuclear spins are each weakly interacting systems whose states may be completely 
specified by the populations of the Zeeman levels. The time evolution of the DNP and 
relaxation processes may then be determined by calculating the rates at which the popula- 
tions of the Zeeman levels of both the paramagnetic electron and nuclear systems change 
as a result of interactions with each other (through electron-dipole-nuclear-dipole coupling), 
with radio-frequency fields, and with the phonon bath of the lattice. In solids, however, it 
is found that the assumption of weakly interacting spins within a spin system cannot always 
be made and therefore that conclusions concerning magnetic processes in solids resulting 
from a rate -equation approach may be erroneous. 

The spin-temperature theory includes the effect of the strongly interacting spins as 
well as the Zeeman interactions and the lattice in a relatively simple way. Within certain 
limits, to be discussed herein, the Zeeman and spin-spin interactions and the lattice can 
each be treated as a thermodynamic reservoir. Each reservoir will have associated with 
it such thermodynamic quantities as temperature and heat capacity. The reservoirs 
exchange energy with each other through the applied radio-frequency fields, certain dipolar 
interactions, and spin-lattice interactions. The time evolution of the DNP and relaxation 
processes are determined by using the density -matrix approach and calculating the rate at 
which the electron and nuclear magnetizations change as the reservoirs exchange energy. 

A detailed derivation of the equations governing DNP and nuclear spin- lattice relaxa- 
tion by using the spin- temperature theory is made. The results of nuclear spin diffusion 
theory in the rapid diffusion limit as reported in the literature are incorporated, and vari- 
ous important limits of DNP, including a new effect, are discussed. A detailed derivation 
is considered necessary since (1) it is not available in the literature, (2) published equa- 
tions governing DNP and nuclear spin- lattice relaxation are either incomplete or incor- 
rect, and (3) all of the assumptions made in previous works are not clearly stated. 

SYMBOLS 


a average distance between nuclear spins 

b pseudopotential radius, 0.68^— 
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b' pseudopotential radius, 0.68 


\ 1/4 


3 , t .. .. \2 T e 


Cm M^nY 


1 + 

n A e 


2 T 2 e 


c ' - ro(*v.) 77 a! 


n 2 e 


Cj(j = 1,2,3,. . .) = 77 Tr 


[. . .[>, vfti’ 1 ) 


,V(t) 


CjQ = 1,2,3,. . .) = tj Tr 


p r , v( tH) 


,v(t) 


CO n 

relative- heat- capacity ratio, — — 

co 2 

L 


relative- heat- capacity ratio, 


relative- heat- capacity ratio, 


(Pn* a ) ! 




1 a 2 

diffusion constant, approximately -±- — — 

10 T 2 n 


d diffusion barrier radius, a( — 

\ y n 


1/3 


a kZ 


distance between the kth electron and the Zth nuclear spin 


E(t) 


time- dependent dynamic nuclear polarization enhancement factor, 


E 


SS 


steady state nuclear enhancement factor 



leakage factor, 
f e leakage factor, 



g(u>) electron spin- resonance line- shape function 

Ne oo 

Tr I I ^S zj §|(T)e-^r 

1=1 

2n Trsf 

y 

g e the electron g- factor, 

91 Hamiltonian operator 

91' Hamiltonian operator without spin- lattice interactions 

unperturbed Hamiltonian operator 
91** Hamiltonian operator in the rotating frame, R9TR - * 

91 jj nuclear dipole- dipole Hamiltonian operator 

®IL nuclear spin- lattice Hamiltonian operator 

91 jg electron-nuclear dipole-dipole Hamiltonian operator 

91 j 2 nuclear Zeeman Hamiltonian operator 

91 r f radio- frequency- interaction Hamiltonian operator 

91 sl electron spin- lattice Hamiltonian operator 



N n T d 

N e T n 

^n^*e 

N e T n 
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electron dipole- dipole Hamiltonian operator 


»ss 

‘tj^SZ electron Zeeman Hamiltonian operator 

secular electron-nuclear dipole-dipole Hamiltonian operator 
91 gg secular electron dipole-dipole Hamiltonian operator 

9^(m=±l,±2) nonsecular electron dipole- dipole Hamiltonian operator 


N e 

qjOB 2 \ 

^SS ~ Y e L 
k>Z = 1 


-3 R 

r kl B kl 


Hq externally applied direct current magnetic field 

H^ amplitude of rotating radio-frequency magnetic field 

fi Planck's constant divided by 2ir 

I x ,Iy,I z x- , y- , z- components of the nuclear spin operator, ^I r j (r=x,y,z) 


x j - ! x j ± %j 


‘rj 


= \/=T 


rth component of the spin operator for the jth nuclei (r=x,y,z) 


k 

m 2 

N e 

N n 


Boltzmann's constant 

second moment of the electron spin resonance line 
number of paramagnetic ions per unit volume 

number of nuclei per unit volume having spin one- half and gyromagnetic 
ratio y n 


P 


parameter describing the numerator of Egg, eq. (99b) 



P’ 


parameter describing the denominator of Egg, eq. (99c) 


P" = P - P' 

Q parameter defined in the expression for r DNp , eq. (100b) 

Q’ = P' - Q 

1/3 


R 


».2a 

T e 


radius of the shell of influence 


’ (f N e)' 


R(t) 


\ l 


unitary operator, 


distance between the kth and Zth nuclear spin 


r ij > r k l distance between the ith and jth or kth and Zth electron 

s x> s y> s z x-, y-, z-components of the electron spin operator, ^S r j (r 


Sj= - S x j ± iSyj 


rj 

S d 

S° 

S± 

s o 


rth component of the spin operator for the jth electron 
saturation parameter, W d T n 
saturation parameter, W^T e 
saturation parameter, W*T 


value of saturation parameter at A = -o> n in figure 3 


S± = W ± T m 
m A n 


S z B i(r) = . 


, O( 0B . ,OB 


SS 


®zi e 


SS 


= x,y,z) 
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S 2 = S 2 + S 2 + S 2 

T d spin- lattice relaxation time of electron dipole- dipole reservoir 

T e spin-lattice relaxation time of electron Zeeman reservoir 

T n spin-lattice relaxation time of nuclear Zeeman reservoir 

T 2e electron spin- spin relaxation time 

mA _ rp 

1 2e 2 A 2e 
T 2e - 3 T 2e 

T 2n nuclear spin- spin relaxation time 

T™ measured electron dipole-dipole reservoir spin-lattice relaxation time 

d 

T™ measured nuclear spin- lattice relaxation time 

nuclear leakage relaxation time 

nuclear relaxation time by the paramagnetic impurity 

T ez electron Zeeman spin temperature 

t time; primes with t indicate integrating parameters 

qjC). 

U(t) unitary operator, e ^ 

V perturbing Hamiltonian operator 

v lattice volume 

\yO electron spin resonance transition probability per unit time induced by an 

interaction with the applied radio- frequency field 
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w d nuclear spin transition probability per unit time induced by the fluctuating 

z- component of the electron spins 

W* solid effect transition probabilities per unit time induced by the applied radio- 

frequency field 

x,y,z coordinates of spatial axes 

a inverse temperature of the electron Zeeman reservoir 

/3 inverse temperature of the nuclear Zeeman reservoir 

/3g Bohr magneton 

inverse temperature of the lattice reservoir 
y inverse temperature of the electron dipole-dipole reservoir 

y e electron gyromagnetic ratio, assumed isotropic 

y n nuclear gyromagnetic ratio 

A = co e - CD 

€ ^ in the shell-of-influence of the kth electron where the electron summa- 

tion index k is suppressed 

3 2 -i 3 

€ W = -2 Y e y n n sin °id cos e d K 

- 3 J 

r j normalization factor, Tr p 

polar angle coordinate associated with the vector r^ , the z-axis parallel to 


8 



p 

p 1 

pR 

T DNP 

^kZ 


0) 

w e 


density matrix 

- 1 R 

density matrix in the interaction representation, U p U 
density matrix in the rotating frame, R(t)p(t)R~l(t) 
dynamic nuclear polarization pump time 

azimuthal angle coordinate associated with the vector rj^ , the z-axis parallel 
to H 0 

angular frequency of applied radio- frequency field 
electron Larmor frequency, - y e HQ 


co L = 


Tr (^ss) 

h 2 Tr S 2 


1/2 




nuclear Larmor frequency, -y n HQ 


W 1 = * y e H l 

Mathematical conventions: 


[a,b] indicates the commutator of the two operators A and B, AB - BA 

An arrow — over a symbol indicates a vector quantity 


Tr 


indicates the trace or sum of diagonal elements 


A tilde ~ over an operator A(t) indicates that A(t) = e 


-iW®t 


SS 


Ae 


SS 


<A> 


indicates a canonical ensemble average of the operator, A, Tr pA 
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BASIC ASSUMPTIONS AND MODEL 


Consider a paramagnetic system consisting of N e v paramagnetic ions (hereinafter 
to be called "electrons") each of gyromagnetic ratio y e , and N n v nuclei having gyro- 
magnetic ratio y n , uniformly distributed in a dielectric, diamagnetic lattice of volume v 
(assumed to be a unit volume here). The Hamiltonian operator for this paramagnetic sys- 
tem, when placed in a homogeneous magnetic field which has a steady component Hq and 
a high-frequency component Hj(t), is given by 

^ = ®sz + ®IZ + ^SS + + ®IS + ®SL + ®IL + ^rf ^ 

where the terms are the electron Zeeman, the nuclear Zeeman, the electron dipole- dipole, 
the nuclear dipole- dipole, the electron- nuclear dipole- dipole, the electron spin-lattice, the 
nuclear spin-lattice, and the radio-frequency interaction Hamiltonians. Explicitly, they 
are: 


N e 

®SZ = X V 

]=1 

N 


( 2 ) 


n 




IZ 


- 2 V 5 0 


(3) 


Z = 1 






N e 






Wss = ry £ r ij 3 

5 • Sj 

- 3r u 2 

(s, • r t j)(Sj • r,j) 


(4) 

i>j = l 






N n 






*H = ?y 1 R kZ 

V T Z 

' 3R kZ 

( T k ’ «b) ( h • 5 k) 

(5) 

s, 

II ' 

I-* 






N e N n 





-1 


®IS y e^nR ^ ^ ^kZ ®k * *Z " 3c *kZ (^k ’ ^kz) ( *Z ' ^kz) 

1-17 1 ' — * 


k=l Z=1 


( 6 ) 
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( 7 ) 


N e N n 

9J rf = -r e R I ^(t) ♦ S k - y n fi ^ H x (t) -t z 
k=l Z = 1 

where sj(]0 is the spin operator of the jth electron (Zth nucleus); y n (y e ) is the nuclear 
(electron) gyromagnetic ratio; K is Planck's constant divided by 2ir; ry is the vector 
distance between the ith and jth electron; is the vector distance between the kth and 

Zth nuclear spin; dj^ is the vector distance between the kth electron and the Zth nuclear 
spin; and H^t) is the high-frequency magnetic field interacting with the electrons and 
nuclei. The 9ig L and terms are discussed in references 10 and 24. The calcu- 

lations to follow will be for electrons having effective spin i and nuclei having a spin i. 

A complete discussion of the concept of effective spin and effective spin Hamiltonians can 
be found in reference 25. 

The following assumptions are now introduced: 

(1) The magnetic field Hq is much larger than the magnetic field experienced by 
each dipole as a result of the effects of the neighboring dipoles. The latter field is called 
the local magnetic field. 

(2) The magnitude of the electron gyromagnetic ratio y e is much larger than that 
of the nuclear gyromagnetic ratio y n , and the inequality N n » N e holds. 

(3) The effective electronic gyromagnetic ratio y e is isotropic. 

(4) The electron and nuclear spin- lattice relaxation times, T e and T n , respec- 
tively, are very long compared with the electron and nuclear spin-spin relaxation times, 

T 2 e and T 2 n , respectively, and there exist times t such that 

T 2e- T 2n <<t<<T e. T n. T d.(wO)- 1 ,(Haw d ) (8) 

where T d is the spin-lattice relaxation time of the electron dipole-dipole reservoir, W° 
is the electron spin resonance transition probability per unit time induced by the interac- 
tion with the applied radio- frequency field, W d is the nuclear spin transition probability 
per unit time induced by the fluctuating z-component of the electron spins, and W* 1 are 
the solid effect transition probabilities per unit time induced by the applied radio- frequency 
field. These terms are discussed in subsequent sections. 

(5) The effects of the lattice may be ignored for times t satisfying the inequality (8), 
and the rate of energy exchange between the various spin systems and the high-frequency 
magnetic field H 1 (t) calculated. The effect of the lattice is then reintroduced. 
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(6) The electron and nuclear resonance lines are homogeneously broadened (refs. 16, 
26 to 28). 

(7) There is no phonon bottleneck (refs. 10, 13, 16, 19, 29, 30 and 31), i.e., the lattice 
heat capacity is assumed to be infinite. 

(8) Additional assumptions will be introduced where appropriate. 

Consider the large constant magnetic field Hq to be applied along the z-axis in the 
laboratory reference frame. Then 


N e 

^SZ = "^e^O X Sz i 
1 = 1 

and 


Define 


^n 

®IZ = _ yn RH 0 ^ l z l 


1 = 1 


w e ~ "^e^O 

w n = -Vn H o 
N e 

s z = ^ s zj 

i=i 


lz _ ^ lz l 

1 = 1 


( 9 ) 


(10) 


(ID 

( 12 ) 


(13) 


(14) 


where w e and w n are the electron and nuclear Larmor frequencies, and S z and I z 
are the z- components of the total electron and nuclear spins, respectively. Since 
|>'e| >> l> / n|’ then |^e| >> l u, nh thus > oscillating field Hj(t) has a frequency w of the 
order of a> e , then the direct interaction of the nuclear spins with this field is negligible. 


12 



This oscillating' field is applied in a plane perpendicular to Hq, so that equation (7) 
becomes 

SW r f = ") / e RH l( s x cos + s y sin wt ) 

where S x and Sy are defined in a fashion similar to S z in equation (13). Defining 
u>i = -ygH^ and noting the identity 


-icoS z t +ia)S z t 
e S x e = S x 


cos cut + Sy 


sin cot 


yields 


t -icoS z t +icoS z t 
- Kcoje S x e 


(15) 


Using equations (9) to (15) with equation (1) gives 


<&' = w e S z + co n I z + 9i ss + Wjg + W n + coje S x i 


-icoS z t g e icoS z t 


(16) 


where the prime is used to indicate that the spin-lattice interactions have been omitted in 
accordance with assumption (5). Note that the Hamiltonian has been written in frequency 
units, that is, fi = 1. 

The time evolution of the system described by the Hamiltonian can be determined 
from the equation of motion of the density matrix p for the system. The density matrix 
obeys the relation (ref. 32) 

3T - < 17) 

again noting that K = 1. 

The density matrix may be transformed by the unitary operator R(t) to become 

p R (t) = R(t)p(t)R -1 (t) (18) 

where p R (t) is the density matrix in the rotating frame and 


R(t) = e icoSzt 


(19) 
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This transformation corresponds to a transformation to a frame of reference rotating 
about the laboratory z-axis with frequency u>. Discussions of transformations to rotating 
frames can be found in references 20, 28, 33, and 34. 

The equation of motion for the density matrix then becomes 



where 


91^(t) — AS Z + w n I z + cOjSjj + R(t)^9Jgg + 91jg + 9Jjj^R ^(0 


and 


A = o> e - a> (21) 

The operator is not altered by the unitary transformation since 



The operators 9igg and 91 jg are of the form 


R(t)9iggR *(t) - '91 gg + ^ 91 ss (23) 

m 


R(t)W^R _1 (t) = 91§ + J (24) 

m 

where the sums are over the values m = ±1, ±2. The operators 9fgg and 91^ are 
the secular parts of the dipolar Hamiltonians, that is, those parts which commute with 
R(t): 



= 0 


(25) 


m m 

and 91gg and 91 jg are the nonsecular parts of the dipolar interactions. The nonsecu- 
lar terms in equations (23) and (24) oscillate at frequencies ±o>, ±2o> in the rotating 
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reference frame and may be ignored. A complete discussion of the effects of these terms 
can be found in references 20, 28, and 33. Explicitly, the secular dipolar interactions can 
be found by carrying out the operation indicated by equation (25) with the result that, in the 
notation of appendix A, 


and 


where 


N £ 


Jigs 


j r kZ ( A kZ + B kz) 


k>Z = 1 


(26) 


N e N n 


^IS ^ X ^kZ^zk^Z + ^kZ^zk^ + e kZ^zk I z) 


k= 1 Z = 1 


(27) 


e kZ = VeVn d kZ (* ' 3 cos2(9 kz) 


(28) 


e kZ = -fve^n sin ^kZ cos ^kZ e ^ d kZ* 


(29) 


and the star indicates complex conjugate. The Hamiltonian in the rotating frame is 
reduced to 


W = AS Z + a> n I z + + Stygg + (30) 

Now, define 

O 

9>o = ^S z + co n I z + 9Jgg (31) 


and 


V = 


ajjSx + 9» IS + 9fjj 


(32) 


so that 
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(33) 


9f R = 3J 0 + v 

It is assumed that (^o) >:> (v) so that V may be treated as a perturbation. Trans- 
forming to an interaction representation by the unitary transformation 

p*(t) = U _1 (t) p R (t) U(t) (34) 

where 

-i<3Ut 

U(t) = e 0 (35) 

results in 

^ = i[p ! , (36) 

where 

V X (t) = U -1 (t)VU(t) (37) 

Equation (36) may be integrated and then iterated to obtain the form 

p'ft) = p 1 © + 2 «>" J ‘ *1 X tl *2 • • • l'"' 1 dl n 

n= 1 

x[[.. p(0), V^tn)], V^tn.ijJ , . . . , v\t{) (38) 

and this result may be differentiated to obtain 
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The basis of the spin- temperature theory is to assume at this point that the density 
matrix in the rotating frame will have the canonical form for times t satisfying equa- 
tion (8). This means that, for the assumptions made thus far, a temperature can be 
assigned to each of the constants of the motion, and these temperatures will represent the 
canonical distribution achieved by each of the commuting terms in * or the times t 

of equation (8). Thus 

expfa(t)AS z - j3(t)u> n I z - y(t)9&] 

p R (t) = — L- j-Y (40) 

Tr / expUa(t)AS z - /3(t)u) n Iz - | 

The inverse temperature a is defined as 


a = 


1 

kT ez 


where k is the Boltzmann constant and T ez is the electron Zeeman spin temperature. 
The inverse temperatures /3 and y are similarly defined. 

The physical model utilized by the spin- temperature theory of DNP and nuclear spin- 
lattice relaxation is that the electron Zeeman, nuclear Zeeman, and secular electron 
dipole-dipole systems in the rotating frame, corresponding to the terms AS Z , a> n I z , and 

in equation (40), will realize a canonical distribution in a time of the order of the 
electron and nuclear spin- spin relaxation times. Therefore, these systems can be treated 
as thermodynamic reservoirs having thermodynamic parameters such as temperature and 
heat capacity. These reservoirs exchange energy with each other through the perturbing 
terms of V. This mixing is realized if the operators represented by the reservoirs do 
not commute with perturbing terms. For times t indicated in equation (8), the rate of 
exchange of energy between the various spin reservoirs is calculated by using the density 
matrix approach. This will result in a time evolution of the temperatures of the spin res- 
ervoirs. The effect of the lattice reservoir on the evolution of the spin temperature is then 
included. A pictorial representation of the reservoirs considered and the various energy 
exchange processes is given in figure 1. The heat capacity of the lattice reservoir is 
assumed to be infinite in comparison with the heat capacities of the spin reservoirs. When 
this situation is not realized, a phonon bottleneck may occur. The effect of the phonon bot- 
tleneck phenomenon on the DNP and nuclear relaxation processes is treated elsewhere 
(refs. 10, 13, 16, 19, 29, 30, and 31). 


17 




Figure 1.- Schematic representation of the thermodynamic reservoirs and 
interactions governing dynamic nuclear polarization and nuclear spin- 
lattice relaxation. 


The electron- nuclear dipole- dipole reservoir and the nuclear dipole- dipole reser- 
voir, as well as the off-diagonal part of the density matrix, have been ignored in equa- 
tion (40). The electron- nuclear dipole-dipole reservoir (that is, the S z ]jl z ^ term in 
equation (27)) contributes to the inhomogeneous broadening of the electron spin- resonance 
line and is assumed to be negligible (ref. 26). This term also gives rise to a barrier to 
spin diffusion, and this barrier is treated phenomenologically when spin diffusion is taken 
into account. The off-diagonal terms in 91^ are retained in V as a perturba- 

tion. The nuclear dipole- dipole interaction may be ignored since its energy (typically 
~5 x 10-29 j) i S V ery small in comparison with those of the electron Zeeman, nuclear 

Zeeman, and electron dipole- dipole interactions in the rotating frame (typically 
-27 

~5 x 10 J). This term is important only when spin diffusion is taken into account. 

The perturbation V(t) reduces to 


N Q N. 


e n 

W l (t) = uijSxft) + w°(t) = Wl Sx(t) + y y (i 

1 _ 1 1 •* ' 


e kZ S zk^ 


+ e kZ s zk5 


k=l 1=1 
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The size of the contribution made by the off-diagonal part of the density matrix is 
discussed by Goldman (ref. 20), and for the time scale considered here, the off-diagonal 
part may be ignored. The validity of the assumption of a canonical form has been shown 
(ref. 35) to be correct within the conditions assumed here. Experimental evidence has 
been published (refs. 20, 33, 36 to 40) verifying the assumption of a canonical form for the 
density matrix for nuclear and electron spin systems. With these approximations, 


p x (t) = p R (t) 

The representation used for the electron system is one in which 31gg, S^, and S z 

are diagonal. Since the nuclear dipole- dipole interaction has been neglected, the nuclei 

are considered independent particles, and the representation used is one in which all the 

o 

single particle operators Ij and I are diagonal. 

The high temperature approximation is assumed to hold for this calculation, that is, 


<a(t)AS z > , 


<0(t)co n I z >| , | <y(t)9fgg> | « 1 


Then, 


P R (t) = 77" 1 


S 

' 1 


1 - «(t)AS z 


- /3(t)o) n I z - y(t)3i 



(41) 


where r\ is the normalizing factor. 

Note that the electron and nuclear magnetizations along the z-axis are proportional 

to the ensemble averages of the operators S, and I , respectively. Using the density 

Z z 

matrix, these values are found as 
<S Z > = Tr p R (t)S z 

and 

<I Z > = Tr p R (t)I z 

where ( ) indicates the canonical ensemble average value and Tr indicates a sum 
over the diagonal elements in a space that includes both the electron and nuclear eigen- 
states. It follows that the time rate of change of these expectation values are (in the high 
temperature limit) 
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(42) 


d <s z > 

dt 


Tr^-S 
Tr dt bz 



Tr S 


2 da 
z dt 


and 


d Pz ) 

dt 


Tr 


dp R (t) 

dt 


I z - 


Trl^ 

z dt 


(43) 


Equations (42) and (43) indicate that the changes in the magnetization are propor- 
tional to the changes of the corresponding inverse temperatures of the reservoirs. The 
calculation of the time evolution of the inverse spin temperatures are derived in the next 
section. 


TIME EVOLUTION OF INVERSE SPIN TEMPERATURES 


Derivation of Equation of Motion for Electron Zeeman 
Inverse Spin Temperature a(t) 

The time evolution of the electron Zeeman inverse temperature a(t) is determined 
by combining equations (39) and (42); that is, 


- A Tr 


si 


) da 
/ dt 


i Tj Tr[p R , V(t)js z 




Sz 



where the superscript I, indicating the interaction representation, has been dropped and 
pR = pR(0). (Henceforth, a time- dependent operator will be explicitly stated.) On the 
assumption that the trace and the time integrals commute, the following terms must be 
evaluated: 
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Using the identity 

Tr [A, B]C = Tr A[B, c| 

yields 

C x = V Tr[s z , p R ]v(t) = 0 

since 



Using the same identity allows C 2 to be rewritten as 

C 2 = V Tr {J R , V(t')j j ^ V(t), S z J 

The density matrix p R (t) can be written in the following way 

P R = i? -1 jl - (a - y) A S Z - (/3 - y) a) n I z - y-3io| 

Inasmuch as 

V(t) = cojSxlt) + 3f^(t) 


(45) 


(46) 


(47) 


(48) 


(49) 
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where S Xj y(t) and ®^(t) are in the interaction representation, C 2 becomes 
C 2 = -o)jA (a-y) Tr S y (t') S y (t) 

+ i Wl u> n 0 - y) Tr S y (t) [l z ,3^(t’)| 

+ i yc^Tr S y (t) [» 0 , V(t’)] 

The second term is zero as can easily be seen by using equation (46) and noting that 

|Sy(t), Izj= 0 


Using the fact that 

i Tr S y (t) [® 0 , V(t)] = Tr S y (t - t') V 
and integrating C 2 as noted in equation (44) yields 



-a^A(a - y) dt' Tr S y S y (t - t') 
+ ya > 1 dt' ^7 Tr S y (t - f) V 


The second term in equation (51) is also zero. This result is shown by letting t = t 
then 


I„ l dt ’ W Tr { S y (t - t,)V } ■ dT Tr["lSxSj,<T) * S y M«°| 


But 


Tr S y (r)^~ = Tr S y (r) 


N e N n 

l !(■ 

k= 1 1=1 


e kl S zk l l + € kZ S zk I Z ) 


= 0 



since ^ has no diagonal elements. The S zk I zZ terms in 9^ have been previously 
dropped but the result shown in equation (53) would not be changed if they had been 
included. Note the identity 

Sy(r) = -S x (t) sin At + §y(r) cos At (54) 


where 


^x,y( T ) _ e 


iSbPqT 
SS s 


x,y 


-i9J 

e 


O 

SS T 


(55) 


The remaining term in equation (52) then becomes 


-u>i 



[jr S x Sy 




cos At - Tr S x S x (t) sin 



But 


Tr S x Sy(T) = 0 


since a rotation of ir about the y-axis sends S x into -S x and leaves Sy(r) unchanged. 
Since the trace is invariant the result follows. Then 


■“1 r dT Tr [ Tr ^ W ] - “1 1 dT Tr Tr sin At 


The term Tr S x S x (t) is proportional to the correlation function of the transverse mag- 
netization, with a correlation time ~ T 2e< Since t is assumed to be very much larger 
than T 2e , the upper limit of the integration may be taken to infinity; thus, 


‘ W 1 




|Tr S x Sy(T)j 


= [Tr S x §x(t) sin a tJ|^_o = 0 


since the correlation function is zero at t = 00 . Equation (51) is reduced to 



r° 


-wJA(a-y)\ dt’ Tr S y Sy(t - t') 


23 



If the same arguments that were previously used are applied, 

J dt' C 2 = -co^A(or- y) ^ dr Tr SySy(r) cos At (56) 

Equation (56) will be further discussed after C 3 and C 4 are deduced. The parameter 
C 3 may be calculated by using some of the relations developed in finding C 2 . Expanding 
Cg and performing the obvious commutation relations results in 

c 3 = -o>2a(o!- y) Tr Sy(t") [®^(f), Sy(t)| 

+ i« n wj (|8 - y) Tr [l z ,3»g(t")] |s x (t’), S y (t)] 

+ icuj 0 .n 0 - y) Tr |l z ,9^(t")] ^g(t'), S y (t)] 

+ iyo.j Tr [W 0 , V(t")] |v(t’), Sy(t)] + 0(o>3) (57) 

The first term in equation (57) can be written as 

Tr Sy(t M ) [®g(t'), Sy(t)]= Tr [s y (t), S y (t")] Swg(t') 

This term is zero for the same reason that equation (53) was zero. The second term in 
equation (57) can be written as 

Tr [l z ,3^(t")] [S x (t’), S y (t)] = Trjjsxa'), S y (t)j, I z ] 9i°(t") = 0 


since 

| x (t’),S y (t)],I^= 0 

The third term in equation (57) is easily seen to be zero by applying the unitary oper- 
ator e * 77 ^ 2 to the argument of the trace. Since S z commutes with I z and < 3 !^(t), the 
rotation sends Sy(t) into -Sy(t) with the result that 
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Tr I Z} Wg(t") wg(t'), S y (t) = -Tr I z ,W°(t") 9fg(t’), S y (t)~ =0 


(58) 


The fourth term in equation (57) when integrated over t" becomes 

iy Wl y dt” Tr [W 0 , V(t”)] [V(t’), S y (t)] = yc^ Tr (v(t’) - V(0)} [v(t’), Sy(t| 

= ycdj Tr [y(t‘), V(t’)| S y (t) 

-ywi Tr V(t’) [s y (t), v] (59) 

The first term in this equation is obviously zero, and the second term can be expanded to 
give 

Tr V(t') |Sy(t),'Vj = Tr ^ Sy(t), SxJ < 3^jg(t') + |Sx(t')> Sy(t)J 
+ Tr Sy(t) jg(t') + 

The first term on the right-hand side is zero for the same reason that equation (53) was 
zero, and the second term is zero when the previous argument leading to equation (58) is 

Q 

used. Therefore, to order of u^, Cg = 0. 

The parameter C 4 can be evaluated by using arguments identical with those used 
in evaluating Cj, C 2 , and Cg, along with the condition that the terms of the order of 
and and higher are neglected. The result is 



It is now assumed that each nuclear spin interacts with only one electron. For dilute 
systems, this approximation is an excellent one. Thus, the sample is broken up into N e 
equivalent "shells- of- influence," each containing one electron and N n /N e nuclei. In 
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this approximation, 



( 61 ) 


where the summation is over one shell of influence, and corresponds to the jth proton 
referenced to an electron at the center of the shell. Integrating C 4 as indicated in equa- 
tion (44) yields 


N n /N, 


Jc 4 dt = -iu> 2 fa - y) Tr £ |e;| 2 iflj £ dt' £ dt" £ df” 

j = l 


x ®x(^ m ) s x (t) 




Nn/Ne 


’AW ^ |. 

+ 1<V) *(0-y>Tr l |«j| 2 lf^ dfj^ dt”J o df 
j = l 




- e 


(62) 


Equation (62) can be integrated as carried out in appendix B. If it is assumed that | ej | 
may be replaced by an average over the shell- of- influence of the electron, that is, an e 2 
such that 


N n /N e 
_ Ne V 
N n A 


j = l 


(63) 


then, 
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jc 4 dt= -A(a-y)|sJeiw°(A) (tt S 2 ) 

+ Ma - y) ^ jw" (cu n - a) + w + (cu n + A) 
- cu n (|3 - y) ^ |w" (cu n - A) - W + (cu n + A) 


Tr S 


2 

z 


Tr S 2 


where 


W 0 (A) = 7T0J J g(A) 


W ± (co n ± a) = 7r(e 2 /cu 2 )cu 2 g(cu n ± a) 


The electron spin resonance absorption line shape function g(cu) is defined as 


Tr \ dT S x (t) S x cos cut 
Ja 

g(w) = 5 

TTTr S 2 


( 64 ) 


(65) 

( 66 ) 


(67) 


and normalized so that 


£ 


dcu g(cu) = 1 


( 68 ) 


Since Sx commutes with neither the A nor the B terms in £bgg (see eq. (26) and appen- 
dix A), S x (t) and thus g(cu) depend on both the A and B terms, as discussed in refer- 
ence 28. The relation (see appendix C) 

poo ~ poo ^ 

Tr L + l7 \ dT S x (t) S x cos cut Tr \ dT S x (t) S x cos cut 
J J Jq 1 

Tr S 2 2 Tr S 2 
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was used to obtain equation (64) from equation (62). The W® and W* are the transition 
probabilities per unit time associated with the allowed electron spin- resonance transition 
and with the solid- effect transitions, respectively. Terms higher than C 4 are of order 
,3 a)? ,3 0)1 

— q> o’ an ^ higher, and are therefore neglected. 

0)3 0)3 A 3 A 3 


Combining equation (44) with equations (56) and (64) to ( 68 ) leads to 


3F = -M^' W ° (i)< “-) 

n j 


|iW( Un -i)^-A ff .( Wn -A) y j 


" ^ Z W+ ( w n + A ) \?nP + A « - ( u n + A )y] 


(69) 


L N n 2 £ 2 

where ^ — is a correction factor to the allowed transition probability per unit time 
o) n 

due to the C 4 term in the perturbation expansion. It is essentially the wavefunction 
renormalization factor given by Jeffries (ref. 13) and Borghini (ref. 31). Because this 
factor is often negligible in comparison with unity, it is omitted hereinafter. The terms 
a, /3, and y in equation (69) are rigorously a(0), j3(0), and y(0). For times t satis- 
fying equation ( 8 ), they are approximately equal to a(t), /3(t), and y(t), respectively, and 
this approximation is used in subsequent calculations herein. 


Derivation of Equation of Motion for Nuclear Zeeman 
Inverse Spin Temperature /3(t) 

The derivation of the equation of motion for /3(t) is similar to the derivation of the 
equation of motion for a(t). Start with equations (39) and (43): 
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(70) 


Traces analogous to the values of^ C of equation (45) may be defined; thus, 



Both Cj and Cg are easily shown to be zero. Evaluating Cg yields 
Cg = n Tr [p R , V(t' )j [v(t), I z ] 

= ~(/3 -y) u> n Tr j^j;g(t)> 

where the other terms are zero. 
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Using equation (27), carrying out the commutation relations, and integrating C 2 as 
indicated in equation (70), gives 

N e N n 

dt' C 2 = -a> n O - y) Tr ^ ^ € ij e kj *j 

i,k=l j = l 


r 

j -t 


dr S zi S zk (r) e" 1 ^ 7 



The product S zi^ zk (r) is related to the correlation function for the z- component of the 
individual electron spins, which involves a correlation time of the order of T 2e . The 
limits on the integral may then be taken as infinity. This expression reduces to 

J dt’ C 2 = -w n (/3 - y) (rr 1 W d (72) 


where 


^n ^e 


Tr 2 Jj e ij e kj L dT S zi^zk( T ) e n 


w d = l^ j = l i>k=l 


2N n 


Tr S" 


(73) 


In the shell- of- influence model, only i = k terms contribute significantly to equa- 
tion (73), and 


N c 


Tr 


w d _ e£ i=l 
2 


n; 


dr S z i§ Z i(r) e 


■iw n r 


Tr 


where e2 is defined in equation (63). A tempting approximation is 


N e 

V* C°° ~ -io) M r f 00 ^ - iu )„ t 

Tr Zj J oo dj S zi S zi ^ e n - Tr j ^ dT S x S x (t) e 
i=l 

j 

in which case W a becomes (by using equation (67)) 


(74a) 
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(74b) 


W d = we 2 g(w n ) 


The significance of this term for DNP and nuclear relaxation is discussed below. 

The validity of the approximation made in the previous paragraph may be checked by 
recalling that g(u>) depends on both the A and B terms in gtfgg. (See equations (26) and 
(67) and appendix A.) Abragam (ref. 28) shows that the second moment M 2 of a reso- 
nance line is (3/2)2 larger when both the A and B terms are used, as compared with the 
case where only the A terms are involved in the M 2 calculation. If the line width Aco 
is taken to be (M 2 ) 1 / 2 , the A(B) terms are found to contribute the proportions to 

Aw. Now (T 2e )“ 1 is proportional to ao>. Let T^ e and T^ 0 be the contributions to 
T 2e from the A and B terms in ®gg. Then 


Ao> oc 




Abragam 's result shows that 



2 J. 

3 T 


2e 



1 _1 
3 T 


2e 


In equations (73) and (74a), the S z j operators commute with the A terms in £Mgg 
but not with the B terms. Hence, W d can be written as 


W d = ir€ 2 g B (a, n ) 


where 


N P 


Tr 


g 


B ( w n) ~ 


dTS zi sg( T )e" 1WnT 


i=l 


2ir Tr S, 
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i ci< OB 


~B iW ss DT 

Sg(r) = e SS S zi e SS 


N e 

gjOB _ 2 y -3 

^SS " y e r kZ ^ 

k>Z = 1 


and B^ is given in appendix A. Thus, g B (co n ) is related to the correlation time 


-.B 


T 2e . For Lorentzians, (ref. 28) 


g(w) = 


( T 2eA) 

1 + ( T 2e") : 


g B (w) = 


( T 2eA) 


1 + 


(4“) 2 


B 

Since T 2e =3T 2e , g(o>) and g B (u>) are noted to be the same to within a factor *3, 
and thus equation (74b) is correct to within a factor a 3. 

The C 4 is found by using the same arguments used to obtain equation (60); thus 

Ci = -i«jA(«- y) Tr jjs y (t’"), ^(t")],^^)] |®°(t), I z ] 

§y(t ,M ),wg(t ,, )|, S x (f) W°(t), I z ] 


^n ~ y ) Tr 


[lz34 0 s(t ,M )], Sx(t") , Sx(t’) [»°(t), I z ] 


Making the same assumptions which led to equation (61) and integrating as indi- 
cated by equation (70) yields 
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N„/N e 


I dt c ; - ‘ a <“ - ■** I bl 2 ', +i 5 £ *’ £ <■'" £ <«"' 


3 = 1 

x S x (t"’)S (t’)(e la,n(t 


•t) e -iw n (t M -t)j 


- a, n W - y)^ Tr 2 bl 2 £ *’ £ *" £ dt ’" 
3 = 1 

X S y (t")Sy(t')(e ia>n( ^" , ' t> + e -la, n (l "' -t )) 


(75) 


When the integration indicated in appendix B is made and equations (63), (66), and (67) are 
used, equation (75) becomes 


J dt C\ = - A^Tr I*)(W“ - W + )(a - y) + w n (Tr l|)(W" + W + )(/3 - y) 


Combining equations (72) and (76) with equation (70) yields 


f = -W d (/3 - y) - JL W - (co n - A) (w n p - Aa - (c n - A)y) 


(76) 


^ W + (w n + A) (co n £ + Aa - (co n + A)y) 


(77) 


Derivation of the Equation of Motion for the Electron 
Dipole- Dipole Inverse Temperature y(t) 

The derivation of the equation of motion for the electron dipole-dipole inverse tern 
perature y(t) can be accomplished by using the fact that energy is conserved in the 
rotating reference frame for times t satisfying equation (8), provided that 

<?o> » < v > 

where Wq and V are defined in equations (31) and (32), respectively. Then 
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l^-=0 = Tr(iS z + <, n I z + wO)^ 


so that 


: (Tr 4 ) f + « 2 (Tr + Tr (»°) 2 £ = 0 


Tr 1“ N 


Substituting equations (69) and (77) into (78) and using = — £ yields 

Tr S* N e 


+^^-W d (i3- y) + Ai W °(A)( a - y) 


dt N e w 2 
L 


^ ^ W '("n - a )f"nP- a “-K - a )>- 


+ 5^ ^ W+ (“n + A ) ("n 15 + a “' (“n + i )>'} 


where 


2 Tr (^y 

L ‘ B 2 TrS 2 


The term u)^ is related to the line width of the electron spin resonance line. It is found 
to be (refs. 19 and 28) 


“L = 3 M 2 


where M2 is the second moment of the electron spin resonance line. 
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Effects of the Lattice Reservoir 

The effect of the lattice vibrations (phonons) is introduced by assuming that each 
spin-lattice relaxation process takes place exponentially with a characteristic spin- lattice 
relaxation time, and that these processes may be added to the equations of evolution for 
a, j8, and y (refs. 10, 19, 20, 28, and 41). Thus, equations (69), (77), and (79) become 

~ = - W°(A)(« -r)+"W'(w n - A) (w n 0 - Aa - (w n - A)yJ 

‘ X W+ K + A ) ( w n0 + Aa - (u> n + A)yJ - ^ (“ “ %,) ( 80a) 

M = -W d (/3 - y) - W" (w n - A) jw n 0 - AO! - (cu n - A)yJ 

- ^ W + (cu n + a) j w n 0 + Aa - (w n + A)yJ - Y~ (/3 - $ h ) (80b) 


% = W°(A)( a ~ y) + —77 W d (/3 - y) 


dt 




N •»£ 


+ 


N n K 


N f 


w 



W"(« n 



co n /3 - A a - (u> n - 



+ ^ J — w+ K + A ) (jV + Aa " ( w n + A )r] “ (? “ ^l) ( 80c > 

W L 

where T e , T n , and T^ are the spin- lattice relaxation times of the electron Zeeman 
(EZR), nuclear Zeeman (NZR), and electron dipole-dipole reservoirs (EDDR), respectively, 
and is the inverse temperature of the lattice. Note that the electron- Zeeman inverse 
temperature does not approach the inverse lattice temperature, but the "colder" tempera- 
ture ^ 0 L . This effect is a result of the transformation to the rotating frame which 
affects the electron Zeeman interaction, but does not affect the nuclear Zeeman and the 
secular electron dipole- dipole interactions. 
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Nuclear Spin Diffusion 

The nuclear dipole- dipole interaction was ignored in the derivation of equa- 

tions (80). This interaction becomes significant if N e « N n , for there may be large 
regions in the solid where the electron- nuclear dipole-dipole interaction is very small 
compared with the nuclear dipole- dipole interaction. The nuclear magnetization in these 
regions will be spatially transported by energy conserving transitions resulting from the 
terms IjT^ of This process is called nuclear spin diffusion. Only a brief discus- 

sion of spin diffusion in the rapid diffusion limit will be given here. More elaborate dis- 
cussions of spin diffusion are given elsewhere (refs. 2, 20, 28, 42 to 50). 

For dilute samples each electron can be considered to interact directly only with the 
nuclei within its shell- of- influence having an outer radius R defined by 



where N e is the number of paramagnetic ions per unit volume. An important parameter 
is the pseudopotential radius b, which is the distance from the electron at which a nucleus 
has approximately the same probability of being relaxed to the lattice by the electron as of 
being flipped by the terms due to neighboring nuclei, and is given by (ref. 44) 

b = Q-68^) 1 ^ 4 ( 82a ) 


The term C is 


c - To ( v e v n K )‘ 


1 + 


and D is the diffusion constant given approximately by (ref. 45) 


(82b) 


D = 


1 a* 


10 T 


2n 


(83) 


Here a is the average distance between nuclei and T 2 n is the nuclear spin- spin relaxa' 
tion time. This idea is now extended to the W d processes, and a second pseudopotential 
radius b' is defined by 


b' = 



(84a) 
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where 


c ' " j| ('Vn") i 


2e 


1 + “n T i 


(84b) 


A Lorentzian form is thus assumed for equation (74b). Another important parameter in 
spin- diffusion theory is the diffusion barrier radius d, which is defined as the distance 
from the paramagnetic ion at which the magnetic field at the site of the nuclear spin due to 
the ion is equal to the nuclear line width (ref. 20), and given very approximately as 


d = 



(85) 


Nuclear spin diffusion is assumed negligible for r £ d, so that d is the inner radius 
of the shell- of- influence. Finally, the diffusion length (DT n )^ 2 , which is the distance 
spin energy diffuses in a time t = T n , is used. 

The regime of "rapid diffusion" is defined when the following relation holds (refs. 44 
and 47): 


a, b, b’ « d « R « (DT n ) 1//2 (86) 

The physical picture is that within the region d ^ r g R, the relaxation and radio-frequency 
transitions tend to drive the nuclear Zeeman inverse temperature to a steady- state value, 
and at the same time the spin diffusion mechanism attempts to keep the total nuclear mag- 
netization spatially homogeneous. 

In the rapid diffusion limit and the shell- of- influence approximation, e^, defined in 
equation (63), is an angular and radial average of |e j | ^ over the shell- of- influence 
(d s r £ R) which ignores the crystal structure. One finds (refs. 10 to 13, 19, 29, 47, 48, 
and 51) 


2 3 Ovy 8 ) 2 

" l0 dV 

so that T n is given by 



(87) 


(88a) 
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where the first term ^1 /T^ ) is due to nuclear relaxation by the desired paramagnetic 
impurity and can be expressed as 



e 2 T e 

1 + (" n T e ) 2 


(88b) 


The second term (l /T^j in equation (88a) is nuclear leakage relaxation by other impuri- 
ties and by other modes of spin- lattice interaction. The expressions given by equa- 
tions (87) and (88) are the appropriate expressions to be used in equations (66), (74), and 
all three equations (80), where 0 is to be understood as an average value over the shell- 
of- influence. 


Nuclear Spin- Lattice Relaxation 

The measured nuclear spin-lattice relaxation rate (tJJ 1 )”* has yet to be deter- 
mined. In the absence of any radio frequency fields, equation (80b) reduces to 


f --w d 8.-r)-(^)(/»-(y 


(89a) 


where 


2 

W d T2e (89b) 

1 + K^e) 

Physically, equation (89a) states that energy in the nuclear Zeeman reservoir can be trans- 
ferred to the electron dipole-dipole reservoir through the interaction as well as to 
the lattice reservoir at the rate 1/T n through processes described by equations (88). 
Therefore, the time evolution of y must also be taken into account. From equation (80c) 
for no applied radio frequency fields 



The general solution for /3 exhibits two time constants. There are several limiting cases 
for which the longer time constant can be easily evaluated; this time constant will be the 
measured nuclear spin- lattice relaxation time Tjj 1 . One limit is « (1/T n ), in which 
case 
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( 91 ) 


tID __ 

n 


n 


A second limit is the situation in which (I/T^) » (l/T n ) and ^N n Wny/N e ci>L) >> !• 1° 
this case y comes to quasiequilibrium ^i.e. , ^ ~ 0^ in times t such that 


Td/^j\ wd ) 1 « t « (w d ) \ T n (92) 

\ N n w n ) 

Thus, equation (90) may be solved for y, which is then substituted into equation (89a). The 
latter can be rewritten as 




dt 


-.m 




(93) 


where 



_1 

T 


n 


1 + 


W q 

!Vi w dT c 

N e ^L 


If 



« 1, equation (94) becomes 




On the other hand, if 



W^Td » 1 , the result is 



T n 


N e “L 1 
+ N n u>n Td 


(94) 


(95) 


(96) 
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Combining equations (88a) and (94) yields 




W 


1 + N n Ct, n 

XT 2 

N e w L 


W a T, 


( 97 ) 


The last term on the right side of equations (94) to (97) is a result of the coupling of the 
nuclear spins with the electron dipole- dipole reservoir. (See refs. 52 and 53.) 

It should be noted that is often taken to be ^ T e (refs. 16 to 20). Experimen- 
tal results given in reference 54 have shown that, for the system studied, the measured 
EDDR spin- lattice relaxation time T™ was of the order of T™ and about two orders 
of magnitude longer than T e . Often it is found that T™ = T^; however, for the limit 
(^n^n/^e^L)^^ >:> (l/T^), equations (89a) and (90) show that y comes into quasi- 
equilibrium with ;3 with the shorter time constant approximately equal to 
(( N n w n/ N e w L) wd ) * > and that /3 and y relax together to the lattice with the longer 

time constant T™ » Tjp, where Tjj 3 is given in equation (96). Thus, the electron 
dipole- dipole reservoir must be included in order to account properly for the nuclear 
spin-lattice relaxation rate. 

THEORETICAL DYNAMIC NUCLEAR POLARIZATION RESULTS 


As was assumed previously, the electron Zeeman and electron dipole-dipole reser- 
voirs are assumed to come to quasiequilibrium [(da/dt) « 0 » (dy/dt)] in times t « T™, 
so that equations (80a) and (80c) may be solved for a. and y, and substituted into equa- 
tion (80b), with all terms retained. This solution is valid provided that (1/T e ), 

(1/T(j) » (l/T n ), W° » W*, and N n /N e » 1. After a tedious calculation, the dynamic 
nuclear polarization enhancement E = can be written as 

^=4;[ E(t) - E ss] <98) 

where r^p is called the DNP pump time and the sample average, steady- state enhance- 
ment Egg is 


E 



(99 a) 


where 
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P = 1 + S d (l + f d c n ) + S°(l + c a r) 


+ S‘ 


o>p co n + A 

— + (• f e + f d c + 

o, n w n 6 


+ s 


COp u> n - A 

— + — — + f e + f ri C. 
aj n a) n e d 


[s°(s d + S + + S“) + f e (S + + S")S d ] 


c 0 * o) n A 

5;T ,U “V' 


Ml 


(99b) 


and 


P’ = 1 + S d (l + f d C n ) + S°(l + c a r) 


+ S + (l + f e + f d C+) 


+ S“(l + fg + f d C_) 

+ [s 0 (s d + S + + S - ) + fg(S + + S“)S d ] [c a R + 1 + f d C n ] (99c) 

Use is made of the saturation parameters S d = W d T n , = W^T e , = W ± T n ; the rela- 

o/o o / 2 2/2 

tive heat- capacity ratios C n = cu n /a> L , C A = A /a> L , C± = (u> n ± A) / co L ; the DNP 

leakage factors f e = N n T e/ /N e T n , f d = f e R; and the electron relaxation ratio R = T d /T e . 
Negligible terms are omitted in equations (99b) and (99c) and in equations (100b), 

(101b), and (102b) that follow. The DNP pump time t dnp is 


1 = _L El (100a) 

r DNP T n Q 


where 



(100b) 


Q = 1 + s d (f d c n ) + S °(1 + c a r) 

+ S + (f e + f d c+) + S“(f e + f d c.) 

+ [s°(S d + S + + S“) + f e (S + + S‘)S d ] f d C n 

and P' is given by equation (99c). 

Alternative forms for Egg and t dnp are easily derived from these equations. 
Egg can be rewritten as 

E SS = 1 + fr (101a) 

where 


+ 



+ S 1 


+ S") + f e (S + 


+ 



(101b) 


and P' is given by equation (99c). When w e » a; n , as is usually the case, w in 
equation (101b) can be replaced to excellent approximation by o> e . Note that Egg = 1 
when S^ = 0 = S + = S“ , as expected. The [s - S’*] term in equation (101b) is DNP by 
the traditional solid effect (SE), discussed in references 10 to 14, 19, 20, 29, 31, and 44, 
and in the special Case 1 that follows. The remaining terms in equation (101b) are DNP 
through EDDR processes, and are discussed in a subsequent section. The DNP pump time 
t DNP can k e rewritten as 


1 

t DNP 


T n 



(102a) 


where 


Q' = [s d + S + 



S°(S d 


+ s + + S' 


') + f e (S + + S')S d |c a r + 1 


(102b) 
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reduces to the T™ 


and Q is given by equation (100b). It is easy to see that Tj^p 
expression in equation (94) when = 0 = S + = S . 

Equations (98) to (102) give complete expressions for DNP and are not to be found in 
the previous literature. It should be noted that the nuclear spin-lattice relaxation time 
T n used in S^, S*, f e , and f d , is not necessarily the measured nuclear relaxation 
time T™. The times T n and T™ are the same in the rapid diffusion limit only when 
the effect of the electron dipole-dipole reservoir is negligible, that is, » 0. 

The wave function renormalization factor (2N n e ^/N e a>n) has been dropped from the 
term in equation (69). This factor would also appear in other terms if the calculation 
were to be carried to higher order. This factor is now estimated for the materials yttrium 
ethylsulphate (YES) and lanthanum magnesium nitrate (LMN), which are often doped with 
paramagnetic impurities and used in DNP experiments. (See refs. 10 to 14, 16 to 19, 29, 
30, 31, and 51 for details and the values used below.) Using equations (12), (81), and (87), 
and the relation g e /3g = y e K, where g g is the electron g-factor and /3g is the Bohr 
magneton, one shows that 


2Nnef 

N e“n 


4tt N n 
5 d 3 




(103) 


or 


2N n e 2 

N e “S 




(104) 


where K = 6.06 x 10"^® and 8.56 x 10“^ CGS units for YES and LMN, respectively. 
The smallest possible value is assumed for d, namely the distance from the electron to 
the nearest nucleus fd = 0.31 nm (for YES), 0.436 nm (for LMN), for protonsj. It is then 
easy to calculate ^2N n e 2 /N e u> 2 j for any case of interest, and, in particular, to find that 


2N„e 


n 


XT 2 

^e w n 


< 0.18 


(105) 


for YES and LMN if g e = 2 and Hq > 2.39 x 10^ A/m (3000 Oe). Finally, note that 

2N n e /^e^n = 2f e i n the limit (a> n T e )^ >3> * anc * (fn/^n) <<; 1* (See equations (87) 
and (88).) 
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Now the results just given are discussed for the limit of well- resolved solid- effect 
transitions (o> n » o>l), in which no more than one of the three transition probabilities 
W®, W* may be nonzero for a given value of A. Three cases are considered: 

(1) (c*> n ± A) = 0, (2) A »0, and (3) (u> n ± A) « 0. Other limits, such as the important 
case of unresolved (o> n % solid-effect transitions, have been discussed extensively in 
the literature (refs. 16 to 23) and are mentioned only briefly herein. 

Case 1: w n ± A = 0 

For Case 1, look first at the DNP exactly at the peak of the solid-effect transitions 
W ± ; that is, o> n ± A = 0. Then, W + = 0 = W®. It is useful to write equation (94) in the 
form 


r m 

n 


_ 1 _ 

Tn 


l + s d (i ♦ f d c n ) 

1 + 


For this case, equation (99) becomes 


(106) 


X + S d [l + f d C n ] + s* 

( T s) + te 

+ f e S ± S d 

r 

/»e\ 

\a> n / 

A 

Uy 

| R + 1 + f d C n 

1 + S d [x + f d C n ] 

+ ^[l + f e ] 

+ f e S ± S d 

C n H 

+ 1 + f d C n 


This equation can be rewritten by dividing the numerator and denominator by 

[l + S d (l + f d C n )j and by rearranging terms using equation (106), the definitions 

SfJb = W*T“ = and fg 1 = fe^n/rJJ 1 ), and the relation fgS* = f“S^. Then, 


1 + Sm 


J SS 


* ^ + f m 

! *>n e . 


1 - # - «?] 

For the usual limit (o> e /a> n ) » 1, this equation can be expressed as 


(107) 



( 108 ) 


Similarly, 
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1 1 [i + s£.(» ♦ '?)' 


t dnp t“ 


1 + fi&t? 


DNP T n 


L- 1 + !ffi 

’? 1 + 


These equations for the pure solid effect agree with references 10 to 13, 19, 29, and 31, 
but include the important EDDR contributions to T™ and thus f|° that these references 
did not include. 

In the limit of saturating microwave power (S^j » l), the DNP enhancement Egg 
becomes 


. 5 2s 

E — i _i_ \ — nj 

ss i * f ” 


and the ideal solid-effect enhancement Ej is attained in the limit of negligible DNP leak- 
age ^f™ « l): 

E ss - E 1 = (* S) (112) 


In equation (112), the 1<< ( CtJ e / a, n) is dropped. Equation (112) is true only in the high - 
temperature limit, since, in general (see references 11, 13, and 19), 


_ tanh (Ko) e /2kT) 
* tanh (Kw n /2kT) 


The DNP leakage factor f™ reduces Egg from E^, as in equation (111). Finally, note 
that t dnp approaches T“[f™/(l + f^)J at S^, » 1, and reduces to T™ 

at = 0. 

This limit (w n ± A = 0) of DNP by the solid effect has been discussed theoretically 
(refs. 10 to 14, 19, 20, and 31) and has been observed experimentally (refs. 10 to 12, 29, 
and 51), and is the case for which the results of the rate-equation approach (refs. 10 to 14) 
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coincide with the results of the spin- temperature theory (refs. 15 to 23), except that the 
latter provides an additional nuclear relaxation mechanism through the electron dipole- 
dipole reservoir. 


Case 2: A » 0 

For Case 2, we look at DNP near the main ESR line (A ~ 0), and still assume that 
o) n » Thus, W° # 0 but W + = 0 = W . If it is assumed that (u, e /u, n ) » 1, 
equations (99c) and (101) become 


E ss = 1 + 



1 + S d (l + f d C n ) + S°(l + C A R) + S°S d (C A R + 1 + f d C n ) 


(114) 


j J 

If S a is negligible, Egg approaches unity. On the other hand, when S a » 1, the result 
at strong ESR saturation (S® » 1) is found to be 


E 



u> n A 


A + ( T e/ T d ) + f e w n 


which has peak values 


Egg 311 of 


(115) 


c.peak _ 
^SS “ 



at 



+ f, 


„i| 

; n 


1/2 


(116a) 


(116b) 


If ^T 6/ /T d j c.2 » f e o,2, which may not be easy to satisfy when o> n » o,^, and if 

(V T d) = 2 , th ^ 


ir peak 
^SS 


0>n 

= ± Ej — - » Ej 

2 v/2o, l 


(117) 
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This effect has been discussed previously (refs. 16 to 20) for unresolved (u> n « w L ) and 
partially resolved (a> n » solid-effect transitions, but it appears to be relevant for the 
well-resolved case (o> n » u> L ) when S d » 1. The basic idea is that the S® transitions 
cool EDDR, and that the nuclei (NZR) are then cooled by the S d transitions connecting 
these two reservoirs. 

A set of parametric curves obtained by using equation (114) with R = (Td/T e ) = 1, 
with a Gaussian line shape function for the electron spin resonance (ESR) transition W° 
and with the parameters in table I, is shown in figure 2. The effect of a nonzero value 

TABLE I.- PARAMETRIC VALUES USED IN CALCULATING CURVES 

IN FIGURES 2 AND 3 

[The values given are typical of values found for samarium 
doped into lanthanum magnesium nitrate at X-band.J 


Parameter 

Value 

U) e 

277-9100 MHz 

o> n 

2tt • 37 MHz 

<4 

9 o 

(2t r) -18.8 MHz 2 


for S d on the enhancement E gg can easily be seen in figure 2. It has not previously 
been observed experimentally or predicted theoretically for the limit of well- resolved 
solid- effect transitions. It should be pointed out that strong coupling between the nuclear 
Zeeman reservoir and the electron dipole-dipole reservoir, indicated by the size of the 
term S d , can be very effective in polarizing the nuclei. For example, the curve indicated 
by ^S d , Sq, f^ = (10, 1, 0) has a peak enhancement larger than the ideal enhancement, 
which in this case is 246, and for values of ^S d , Sq, f d ^ = (10, 10, 0), the peak steady-state 
enhancement is over 800 (more than 3 times Ej). The effect of leakage (f<j = f e R) is to 
decrease the steady- state enhancement E gg , but even for large leakage factors, polariza- 
tion should still be observable. 


Case 3: w n ± A x 0 

Case 3 is near the solid-effect transitions (o> n ± A ~ o) with S* i= 0, but S® = 0 = S T , 
a> » CJ T . Exactly at the S* peaks, Case 1 ((w n ± A) = 0), EDDR contributes indirectly to 
the DNP only through nuclear relaxation processes, that is, through T n and f g . However, 
slightly off the S* peaks (w n ± A « 0), EDDR may contribute explicitly to the DNP. Using 
equations (99c) and (101), yields for this case (taking o» = o> e ) 
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Figure 2.- Parametric curves of the steady- state enhancement Egg from equa- 
tion (114), as a function of the applied microwave frequency v = (u>/ 27 t), near 
the main ESR transition W^. The solid-effect transition W”, centered at 
9137 MHz, is well separated from W^. The Egg curves are labelled by 
^S^, Sq, f^). R = = 1. is assumed to be a Gaussian: 

S®(v) = Sq exp Q (i^9100)^/12.5], where Sq is a constant. 
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(118) 




O) 


n 


S* [ + IS ]f l - f e S d Rl 


2 


E SS = 1 + 


Cl) 


1 + S d (l + f d C n ) + ^[l + f e + I d C ± + f/^R * 1 * f d C„) 


If f p S d « 1, EDDR and DNP effects are negligible, as equation (118) reduces to equa- 
tion (108) for all A values near (+ cu n ). However, if f e S a « 1, equation (118) predicts 
asymmetries in Egg as a function of A, as well as slightly higher Egg values than 
those in equation (108). In this case EDDR is cooled by the solid-effect transitions S^, 

j 

and this cooling is conveyed to the nuclei by the S a transitions. This proposed new DNP 
effect involving the terms (fgS^S^ is a higher order analog of the Case 2 (S®)S d mech- 
anism discussed previously for the limit u> n » u^. In the general case of unresolved 
satellite transitions ^co n « both mechanisms operate simultaneously, as is easily seen 

from the Js^(S d + S + + S ) + f e (S + + S )S d J terms in equations (99) and (101); however, 
the well-known S^(S d + S + + S ) mechanism would then ordinarily predominate. The 
two mechanisms are in principle separable for well- resolved satellites, which should be 
useful for experimental tests of them. 


when 

ment 


It is interesting to note that equation (118) reduces to equation (115) with A « + co n 
f e S d » 1 and S* » |jL + S d (l + f,jC n )J, and that this equals the solid effect enhance- 

(+ COp/cOn) _ 

E SS = ^"m — f rom equation (111) exactly at A = •+ u> n . 


Figure 3 presents parametric curves of the steady-state enhancement Egg given 
in equation (118) for various nonzero values of S d and two values of f e = f^ (for the 
ratio R = 1). The electron spin resonance line shape is assumed Gaussian, and relevant 

j 

parameters are given in table I. The curves for S a = O are also shown in figures 3(a) 
and 3(b) for comparison purposes. It can be seen from both figures that even small values 
of S d will distort and broaden the enhancement curve when Sq = 10. The distortion at 
lower levels of saturation is not as evident but there is a noticeable broadening effect 
especially when f^ = 10. Comparing figures 3(a) and 3(b), it is evident that the effect of 
the leakage term is to broaden the enhancement curve for fixed Sq. 


CONCLUSIONS 


A detailed derivation of the equations governing dynamic nuclear polarization (DNP) 
and nuclear spin-lattice relaxation using the spin temperature theory has been carried to 
second order in a perturbation expansion of the density matrix. This theory is valid for 
ionic crystals weakly doped with paramagnetic impurities assuming (1) a completely 
homogeneously broadened electron spin resonance line, (2) an isotropic electron 
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+ 20 ) — 



( a ) fe - f d ~ 


Figure 3.- Parametric curves of E gg from equation (118), as a function of 
v = (co/ 27 t), near the solid-effect transition W~. The main ESR line , 
centered at 9100 MHz, is well separated from W". The E gg curves 
are labelled by ^S d ,SQj. R = (T d /T e ) = 1. W" is assumed to be a 
Gaussian: S~{v) = Sq exp £ (v - 9137 ) 2 /l2.5] where Sq is a constant. 
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+ 4p — 

(1 



9120 9128 9136 9144 9152 

v, MHz 

(b) f e = f d - io. 

Figure 3.- Concluded. 

g-factor, (3) rapid nuclear spin diffusion, (4) the high temperature limit, and (5) no phonon 
bottleneck. Complete expressions for the DNP steady state enhancement E gs , the DNP 
pump time, and the measured nuclear spin- lattice relaxation time, which include effects of 
the electron dipole-dipole reservoir (EDDR) on the nuclei, have been derived. 

Several cases for the limit of well-resolved solid effect transitions have been 
examined. The conventional results for DNP by the solid effect exactly at the solid effect 
transition peaks have been obtained, but with EDDR effects on the nuclear relaxation and 
DNP leakage factor included. It is suggested that direct EDDR contributions to DNP 
excited by saturation of the main ESR line, previously studied for unresolved and 
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partially- resolved solid effect transitions, may be extended to the well resolved case, 
leading to Egg values greater than the conventional solid effect. A new DNP effect is 
proposed: saturation of the solid effect transitions may cool EDDR in favorable cases 

and lead to asymmetries in the shape of the DNP enhancement curve as well as to slightly 
higher peak polarizations than those at the solid- effect transition peaks. Parametric 
curves of these EDDR effects on DNP are given for typical parameter values. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 21, 1974. 
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APPENDIX A 


DIPOLE- DIPOLE HAMILTONIAN OPERATORS 


It is customary to write the magnetic dipole-dipole interaction Hamiltonian for a 
single spin species as foUows: 


N 


®SS - Ve" 2 


2 r k2 f A kZ + B kZ 

7 * 


+ C kZ + D kZ + E kZ + F kZ 


k>Z = 1 




(Al) 


where 


A kZ “ S zk S zZ l 1 " 3 cos2 0 kz) 


BkZ = - 4 (l - 3 cos 2 (sjsf + Sft + ) 


(A2) 

(A3) 


C kZ " - | sin «kZ cos e U e ~ iCP]d { S z^ + S z l*t) 


(A4) 


D kZ ” ” | sin 0 kZ cos 0 kZ e ( S zk S Z + S zZ S k) 


(A5) 


F - 3 Gin 2 0 e~ 2i<p te Q + S + 

E kZ ■ 4 Sin 6 W e ®k S Z 


(A6) 


F kz =-f sin 20 kz e+21 ^ Vz" 


(A7) 


where 6 ^ and are the polar coordinates of the vector r^ , the z-axis is parallel 

to 0 q, and 


®k = s kx ± i& 


ky 


(A8) 


The Hamiltonians gtjg and can be similarly written. 
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APPENDIX B 


EVALUATION OF TRIPLE TIME INTEGRALS 

The triple integrals in equations (62) and (75) are carried out by integration by parts. 
The first term on the right-hand side of equation (62) is integrated to illustrate the method. 
The remaining integrals in equations (62) and (75) are treated in the same manner and lead 
to equations (64) and (76), respectively. 

Define J(o> n ) as 


J(w n ) = f * dt' f 1 dt" e la>n(t _t } ^ dt"' S x (t’”) S x (t) 

Jq v»q Jq 


(Bl) 


and integrate this expression by parts 


i pt f iu> n (t"-t') pt" 
j( w ) = JL\ dt’ ie \ dt ,M S x (t’") S. 

iw n J Q 1 J 0 

- dt” e iWn(tM “ t,) S x (t") S^tjJ 


(t) 


t"=t' 

t M =0 


which reduces to 


J(w n ) = r*- C dt' C* dt" S x (t") S x (t) 

n iw n J 0 J 0 


1 pt pt f ico n (t TT -t T ) 

- rM dt’ \ dt” e nV S x (t") S x (t) 

iu>n j q J 0 


Integrate the second term on the right in the same manner, 


i pt ft' 1 pt iaMt’-t; / v 

J(o> n ) = 4- dt’ dt" S x (t") s x (t) - -M dt’ e n S x (t') S x (t) 

io> n Jq J 0 uit 0 

n 

+ — s (* dt' Sx(t') s x (t) 


(B2) 


(B3) 


(B4) 
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APPENDIX B - Concluded 


The first term on the right of equation (62) is 


F = 


N n /N e 

-Au> 2 (a-y)Tr ) 

j = l 


2 ft 


^J(ai n ) + J(-ai n ) 


2Acu 


<J 0 ‘ 


u N n /^ e t 

k<* - y) Tr £ |ej | 2 1+17 J dt’ S^t'-OS* cos A(f-t) 

j = l ° 


N n^e 

AcJ l V i 

+ — («-y) Tr 2 b 

J = 1 


2 I + L 
J 1 


x 



' Sx(t'-t) Sx |cos (a> n + A)(t'-t) + cos (cu n - A)(t'-tj] 


(B5) 


If equations (63) and (65) through (68) are used, equation (B5) reduces immediately to 
equation (64). 
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APPENDIX C 


TRACE CALCULATIONS 

This appendix illustrates some of the trace calculations used in the shell-of-influence 
approximation. The numerical values cited are for nuclei and paramagnetic ions (called 

electrons herein) of spin—, that is, for I. = ^ and S = ^ respectively. Since there are 

^ J lx it 

two spin spaces involved, the trace must be taken over the product space. For example, 

Tr 1 = Trg • Trj 1 = Trg 1 • Trj 1 (Ci) 

where Trg(Trj) signify the traces in the subspace of the electron S (or nuclei I), and 
where 

Tr g l = 2 • S + 1 (C2) 

and 

N n /N e 

Tr x l = (2*1 + 1) (C3) 


Here the number of nuclei in the subspace of I is N n /N e , that is, the number in each 
shell-of-influence. There is one electron in each shell-of-influence. Other trace iden- 
tities used are 


Tr s S^ = (1/3) S(S + 1)(2S + 1) 

(C4) 

Trjl^ = (l/3)^N n /N e )l(I + 1)(2I + l) Nn ^ e 

(C5) 


and 


Tr^. = (1/3) 1(1 + 1)(2I + l) Nn/Ne (C6) 

These identities are used to evaluate equations involving traces over the product space of 
the S- and the I-spins, for example, in evaluating the quantity Tr L + L F(S)/Tr S^ , where 
F(S) is some function of the S- spins only. This evaluation is carried out as follows: 
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APPENDIX C - Concluded 


TvTtf 


F(S) 


Tr S: 


Trjlfl] - Tr s F(S) 

x 


Trjl 


Tr S S z 


1 

2 


Tr g F(S) 

Tr s S^ 


x Tr g F(S) Trjl 

2 Tr s s z Tr r 1 2 


1 Tr F(S) 

2 Tr 
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